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Abstract 

Several dynamical stages of the Marangoni convection of an evaporating sessile drop are obtained. We jointly 
take into account the hydrodynamics of an evaporating sessile drop, effects of the thermal conduction in the 
drop and the diffusion of vapor in air. The stages are characterized by different number of vortices in the 
drop and the spatial location of vortices. During the early stage the array of vortices arises near a surface 
of the drop and induces a non-monotonic spatial distribution of the temperature over the drop surface. The 
number of near-surface vortices in the drop is controlled by the Marangoni cell size, which is calculated 
similar to that given by Pearson for flat fluid layers. The number of vortices quickly decreases with time, 
resulting in three bulk vortices in the intermediate stage. The vortex structure finally evolves into the single 
convection vortex in the drop, existing during about 1/2 of the evaporation time. 



1. Introduction 

Evaporation of a drop in an ambient gas was con- 
sidered since Maxwell time mainly as diffusion of 
vapor from a near-surface layer [J 0, Q • Classical 
quasistationary theory of evaporation of a drop does 
not include effects of fluid dynamics in the drop and 
only partially takes into account basic effects of heat 
transfer. Recently a number of new applications has 
stimulated much attention to the problem again. 
This is associated, for example, with preparing 
ultra-clean surfaces 0, S 0- 0] , protein crystallogra- 
phy 0, 0] , the studies of DNA stretching behavior 



Other example, where depinning of the contact 
line is revealed, is the effect of evaporative con- 
tact line deposition th e so-called coffee-ring ef- 



and DNA mapping methods ll(X 111. 1211 . developing 
methods for jet ink printing Y&_ 14 Hl5| . and with 
other fields (see, for example. |l6l|). Of particular 
interest is the process of evaporation of the drop 
containing the colloidal suspension [It], 18, 19[ . One 
of typical experimental situations is an observation 
of an evaporating drop lying on a substrate, where 
the contact line of the drop is pinned to the edge of 
the substrate. Ordered structures of nanoparticles 
can arise on the drop surface during the evaporation 
process. After the drop dries out the ordered struc- 
tures are left also on the substrate. An important 
example is the self-assembly of long-r ang e-ordered 
nanocrystal superlattice monolayer [lT], EM H^ |- 
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feet [Io|,00,00,|25|. These aspects altogether 
resulted in the interest to the problem and respec- 
tive significant activity of experimentalists and the- 
orists during the past decade. 

For theoretical description of problems men- 
tioned above, the classical theory of drop evapora- 
tion, which was used up to recent years, turned out 
to be insufficient, too simplificated and not includ- 
ing a number of important processes. For example, 
role of pinned contact line during the evaporation 
process, and fluid dynamics effects within the evap- 
orating drops of capillary size, were analysed and 
new important aspects were found. As a result, a 
development in experimental studies and progress 
in understanding of evaporation process have been 
achieved during past decade [U El 

It was found, in particular, that the evaporating 
flux density is inhomogeneous along the surface and 
has an integrated singularity on approach to the 
pinned contact line (20 . [2lj . The resulting inhomo- 
geneous mass flow can strongly modify the temper- 
ature distribution over the drop surface and, hence, 
the Marangoni forces associated with the tempera- 
ture dependent surface tension. The convection in- 
side a droplet 0, Ez , 
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appears to be quite different from the 
classical Marangoni convection in the systems with 
a simple flat geometry [42j, |43| . Thermal conductiv- 
ity of the substrate can also influence the formation 
of flows within a liquid drop since it is the magni- 
tude of the conductivity which determines the sign 
of the tangential component of the temperature gra- 
dient at the surface close to the contact line and, 
therefore, the direction of the convection [29j |. 
The observation of the distinct stages of the evap- 
11, EF ' 



oration process 



has revealed that the 



longest and dominating regime of the evaporation 
process is the constant contact area mode, where 
the contact line is pinned. On a later stage the con- 
tact line gets depinned and the different regime, the 
constant contact angle mode switches on. Finally, 
the drying mode follows, in which the height, the 
contact area and the contact angle rapidly decrease 
with time. 

Recently, evolved numerical calculations have 
been carried out for providing simultaneous de- 
scriptions of the fluid dynamics, the vapor diffu- 
sion and the spatial temperature distribution in an 
axially symmetrical evaporating sessile drop. The 
cor resp onding quasistationary problem was studied 
in [27L [28} . The dynamical description of the liq- 
uid drop evaporation was developed in [47J . Taking 
into account quantitatively the important compo- 
nents of the evaporation process, we have obtained, 
in particular, the time evolution of the temperature 
and the fluid convection in the drop. 

2. Dynamics of vortex structure 

Our calculations demonstrate the presence of sev- 
eral characteristic stages of the thermocapillary 
convection. During the early regime of the dynam- 
ics of the Marangoni convection, for various liquids 
and drop sizes, the vortices arise near the surface of 
the drop. For the toluene drop, this regime quickly 
arises and evolves up to t « 0.3 s. This is quite a 
short time period for the toluene drop as compared 
with the total evaporation time ~ 550 s, but it ad- 
mits an experimental study. The vortices grow, the 
number of vortices decreases, and eventually they 
evolve into the bulk of the drop. This can be seen 
in Fig. [T] where the vortex structure contains four 
pairs of near-surface vortices and a corner vortex, 
and the temperature displays just four humps at 
the surface. This number of vortices in the drop 
is controlled by the Marangoni cell size, which is 
calculated similar to that given by Pearson for flat 



fluid layers. The existence of near-surface vortices 
and the associated humps in the profile of the sur- 
face temperature, become more pronounced with 
the decrease of the viscosity of the liquid. There 
are no near-surface vortices when the viscosity in- 
creases more than in four times as compared with 
the toluene drop. Near-surface vortices were ob- 
served experimentally in studying evaporating ses- 



31 



sile drops of FC-72 on a copper substrate 3(1 

As seen in Fig. [TJ the extrema of the surface tem- 
perature correspond to the change of sign of the 
tangential component of velocities at the surface. 
The reason for this behavior is that the fluid flow 
moves from the higher to the lower temperature re- 
gions of the surface, because the surface tension de- 
creases with increasing the temperature. The flows 
result in a redistribution of the temperature due to 
the convective heat transfer. The number of the 
surface temperature humps and the number of the 
near-surface vortices decrease during their evolu- 
tion. 

If the thermal conductivity of a substrate is large 
compared to that of the liquid, then the tempera- 
ture can be maintained practically constant at the 
substrate-fluid interface. This is the case, in partic- 
ular, for the silicon nitride substrates used in exper- 
iments [I?], O, Gil ■ The silicon nitride is a material 
with high thermal conductivity, approximately in 
three orders larger than for the toluene. For this 
reason, the boundary condition for the tempera- 
ture distribution at the substrate can be reduced 
to the constant temperature. Heat transfer between 
the substrate and the drop plays an important role 
in establishing the temperature profile in the drop. 
High substrate conductivity also excludes a possi- 
bility for the reversal of the Marangoni flows (29|, 
taking place for substrates with relatively small 
thermal conductivity. 

The initial conditions chosen when generating 
Fig. [T] were: zero velocity values, constant room 
temperature, and absence of vapor. For a descrip- 
tion of a real experiment they have to be slightly 
modified to include the weak stochastic distribu- 
tion of the surface temperature and velocities. We 
have carried out such calculations using a small- 
scale stochastic initial conditions. This has changed 
the particular behavior of the fluid dynamics only 
on the initial stage of the process, where a large 
number of small-scale surface vortices arise. Then 
the surface vortex structure quickly evolves into 
exactly the same one as we obtained for the ba- 
sic initial conditions. This result demonstrates the 
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Figure 1: (Color online.) Left panel: The velocity field at t = 0.16 s. Right panel: Temperature (solid line) and velocitiy 
(dashed line) along the drop surface as functions of the radius at t = 0.16 s. The velocity changes its sign at crossing points 
with the horizontal axis. 
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generic character of the near-surface vortex regime 
in Fig. [T]at the early stage of the formation of the 
Marangoni convection. 

An initial velocity field within the sessile drop can 
also be strongly disturbed right after the drop has 
fallen down on the substrate, or for some other rea- 
son. We model such a situation by choosing random 
initial conditions for the bulk velocity field, which 
are in agreement with the continuity equation for 
the incompressible fluid. We find that strong dis- 
turbances of the bulk velocities can noticeably mod- 
ify the initial stage of the drop dynamics, but they 
do not modify its main stage. For the initial random 
velocities of the order of 5 cm/s (which well exceeds 
the typical velocity in the vortex 1 cm/s), the dif- 
ference between the dynamics of the disturbed and 
resting drops disappears after t w 0.5 s. This ver- 
ifies the stability of the large-scale drop dynamics 
with respect to disturbances of the initial tempera- 
ture and velocity fields. 

During the enlargement of the near-surface vor- 
tices, their number decreases, and the convection 
involves the bulk of the drop. As a result, for 
t w 0.45 s, three bulk vortices control the veloc- 
ity and temperature fields in the drop, as seen in 
left panels in Figs. 12131 During the coexistence of 
three vortices, the corner vortex starts growing at 
the expense of the other two vortices, and eventu- 
ally at t w 2.0 s it occupies the whole drop volume. 
A spatial dependence of the temperature along the 
drop surface is nonmonotonic, if the drop contains 
more than one vortex (see Fig. [3]). Right panels 
in Figs. 12131 demonstrate how in the single- vortex 
regime effects of Marangoni forces drive liquid along 



the surface to the apex, where the fluid penetrates 
along the symmetry axis in the depth of the drop. 

The regime with the single vortex represents one 
of the main stages of the dynamics of the evapo- 
rating sessile drop. It lasts up to t w 250 s. More 
than half of the drop mass evaporates during this 
period of time. If the initial values of mass, height 
and contact angle of the drop are m = 8.7 mg, 
h = 0.1314 cm, 6 = 1.2045, then at the moment 
t = 250 s we find m — 4.0 mg, h = 0.0685 cm, 
6 = 0.716. In particular, h/(2r ) « 0.17, i.e. the 
drop shape is noticeably flattened. We remind that 
the total time of the evaporation of the toluene drop 
is 508 s. 

The quasistationary single-vortex state loses its 
stability at t 250 s and the vortex acquires a pro- 
nounced nonstationary character. During this non- 
stationary regime, the fluid pulsations take place. 
The characteristic frequency of the pulsations cor- 
responds to the circulation period 0.15 s of a fluid 
element in the original vortex. Initially the pulsa- 
tions are concentrated near the center of the orig- 
inal vortex. Then, the single-cell pulsating state 
breaks into two-center (and later three-center) pul- 
sating structure (see Fig. 10 in [47jj ) . Eventually at 
t ~ 300 s a quasistationary state with three vortices 
arises. 

The numerical calculations of the fluid dynamics 
were tested with several different mesh sizes. The 
respective results are qualitatively identical and 
show reliable convergence of the quantitative char- 
acteristics. For example, the single-vortex regime 
was found to arise at 3.48 s, 2.48 s, 2.2 s, 2.06 s, 
2.01 s for 100x100, 150x150, 200x200, 250x250, 
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Figure 2: (Color online.) Left panel: The velocity distribution at t = 0.5 s. The stage of drop dynamics with three vortices 
takes place from t fa 0.45 s to t f» 2.0 s. Right panel: The velocity distribution at t = 30 s. A distribution with a single vortex 
takes place from t fa 2.0 s to t fa 250 s. 




Figure 3: (Color online.) Distributions of temperature within drop for the stages of drop dynamics with three vortices (left) 
and a single vortex (right). The distributions are taken at t = 0.5 s and at t = 30 s correspondingly. Temperature scale shown 
at the right column. 
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300 x 300 mesh elements covering a half of the drop 
cross-section. 

3. Numerical method 

The simultaneous calculation of the physical 
quantities in the drop can be partitioned into sev- 
eral steps: 

1. We apply to the diffusion equation du/dt = 
DAu the implicit finite difference method us- 
ing irregular mesh outside the drop and a vari- 
able time step. We use a boundary interpola- 
tion in a vicinity of the drop surface [47j. For 
the boundary conditions we take u = u s on the 
drop surface, u = far away from the drop, 
du/dr — and du/dz = on the axes r = 
and z = correspondingly. 
In order to calculate accurately the vapor den- 
sity in quite a large region, it is convenient to 
use irregular mesh. We use the mesh with suf- 
ficiently small steps near the drop surface and 
with exponentially increasing steps, which are 



chosen in accordance with the decay of the va- 
por density, which can be estimated analyti- 
cally. 

2. Calculations of the stream function ijj and ve- 
locities v inside the drop are based on Eq. 
d 2 ip/ dr 2 — dip I rdr + d 2 ip/ dz 2 = ry, where 7 is 
vorticity. The implicit finite difference method 
with a regular mesh inside the drop is applied. 
We use a boundary interpolation near the drop 
surface. For the boundary conditions we take 
■ijj = at all boundaries: on the surface of the 
drop and on the axes r = and z = 0. 

3. We solve Eq. dj(r, z)/dt + (v • V)j(r, z) = 
v (Ajir, z) — 7(7-, z) /r 2 ) to obtain the vorticity 
7 inside the drop. The explicit finite difference 
method with a regular mesh inside the drop is 
used. We use a boundary interpolation close to 
the drop surface. For the boundary conditions 
we take 7 = for r = 0; 7 = dv r /dz for z = 0; 
7 = da /{j]ds) + 2v T d(p/ds on the drop surface, 
where da/ds = —a'dT/ds is the derivative of 
the surface tension along the drop surface. 
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For calculating the temperature T inside the 
drop, the explicit finite difference method with 
a regular mesh is applied to the thermal con- 
duction equation dT/dt + v- VT= kAT. We 
use a boundary interpolation in a vicinity of 
the drop surface. The boundary conditions 
take the form dT/dr = for r = 0; T = T Q 
for z = 0; dT/dn = -Q (r)/k = -LJ{r)/k 
on the drop surface. Here Qo(r) is the rate of 
heat loss per unit area of the free surface, n is 
a normal vector to the drop surface. 
During the iterative procedure, the drop shape 
is recalculated in accordance with the evapora- 
tive mass loss for the respective time interval. 
This is possible due to a quasistationary char- 
acter of the shape change of the drop. Devia- 
tions of a drop profile from a spherical cap are 
taken into account [47], 



4. Discussion 

We have developed the approach for studying the 
evaporation and fluid dynamics of a sessile drop of 
a capillary size and applied it for the description of 
the toluene drop evaporation. In particular, the ap- 
proach allows to obtain the time evolution of vortex 
structure and temperature distribution in evaporat- 
ing sessile drops. The appearence of near-surface 
vortices on the early stage of the evaporation pro- 
cess of sessile drops is predicted. We also obtained 
three bulk vortices in the intermediate stage. They 
finally evolve into the single convection vortex in 
the drop, existing during about 1/2 of the evapo- 
ration time. Possible development of the method 
include description of hydrodynamics of colloidal 
drops, and study of influence of substrate conduc- 
tivity on the dynamics of vortex structure. 
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